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动态 选择 策略 的 和 声 教 与 学 混合 算法 * 
黎 延 海 ， 拓 守恒 ， 雍 龙泉 


(陕西 理工 大 学 数学 与 计算 机 科学 学 院 , 陕西 汉中 723001) 


摘 要 : 为 提高 对 多 种 不 同类 型 问题 的 优化 性 能 , 提出 了 一 种 基于 和 上 声 搜 索 和 教 与 学 优化 的 混合 优化 算法 (HHSTL)。 
在 不 同 的 进化 阶段 ，HHSTL 算法 依据 种 群 活跃 率 及 种 群 最 优 个 体 更 新 率 动 态 地 确定 和 上 声 算 法 或 教 与 学 算法 作为 下 一 
周期 种 群 更 新 方式 的 比例 ， 并 在 标准 教 与 学 算法 中 增加 了 “自学 "策略 来 提高 算法 的 全 局 寻 优 能 力 。 通 过 16 个 不 同类 型 
的 Benchmark 函数 进行 仿真 ， 并 与 7 种 优秀 算法 进行 结果 比较 及 Wilcoxon 秩 和 检验 分 析 ， 结 果 表 明 HHSTL 算法 汲取 
了 和 声 搜索 和 教 与 学 优化 算法 的 优点 ， 具 有 求解 精度 高 、 稳 定性 好 等 特点 ， 能 够 求解 更 多 的 较为 复杂 的 优化 问题 。 
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中 图 分 类 号 : TP301.6 — doi: 10.3969/j.issn.1001-3695.2018.06.0399 


Hybrid algorithm based on harmony search and teaching-learning-based optimization using 
dynamic selection strategies 
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c» Abstract: This paper presents a hybrid algorithm based on Harmony Search and Teaching-Learning-based Optimization 
(HHSTL) for solving various types of optimization problems. In different evolution stages, HHSTL algorithm dynamically 
determines the proportion using the harmony search (HS) or Teaching-Learning-based optimization (TLBO) algorithm as the 


means of updating the next generation of population on basis of the population activity rate and population optimal individual 


phe renewal rate, in addition , it proposed a "self-study" strategy to improve the global optimization performance of the TLBO 


algorithm. It employed 16 various types of Benchmark functions to investigate the performance of HHSTL, and compared with 


seven algorithms by Wilcoxon rank sum tests . The experimental results indicate that HHSTL algorithm adopted merits of HS 
and TLBO, and can solve more complex optimization problems because of its high solution precision and excellent stability. 
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在 科学 计算 和 工程 技术 中 ， 经 常会 遇 到 类 型 各 异 的 优化 问 法 在 “ 教 ” 阶 段 存在 固有 的 起 源 偏差 ， 对 以 坐标 原点 为 最 优 解 
Ui, 单一 的 优化 算法 往往 具有 偏好 性 , 适用 的 范围 有 限 , 因此 ， 的 问题 有 着 原始 的 偏好 ， 但 对 最 优 解 不 在 原点 的 问题 则 优化 效 
将 县 有 不 同 优 点 的 算法 进行 融合 能 够 扩大 求解 问题 的 范围 ， 具 果 不 太 理 想 。 这 两 种 算法 具有 不 同 的 特点 ， 适 合 解决 的 问题 类 
有 重要 的 实际 意义 。 和 声 搜 索 算法 (harmony search, HS) DA 型 也 各 不 相同 ， 因 此 将 两 种 算法 进行 融合 ， 能 达到 优势 互补 ， 

是 Geem 等 人 通过 模拟 现实 中 的 乐队 即兴 创作 过 程 而 于 2001 年 各 取 所 长 的 效果 。 目 前 ， 关 于 将 这 两 种 算法 相 融 合 的 研究 主要 
提出 的 一 种 启发 式 优化 算法 。HS 算法 具有 较 强 的 全 局 搜索 能 有 Aydogdu 等 人 "将 HS 策略 融入 TLBO 算法 而 提出 了 混合 算 
力 ， 适 用 于 多 模 态 优化 问题 的 求解 ， 但 局 部 开发 能 力 较 弱 ， 使 法 , 并 应 用 于 钢 架 结构 设计 问题 ; Tuo 等 人 四 提出 了 面向 高 维 复 
得 在 优化 单 模 态 问题 时 所 得 到 的 解 的 精度 不 高 中 。 教 与 学 优化 杂 优 化 问题 的 HSTLBO 算法 , 算法 应 用 了 一 个 自 适应 策略 对 两 
算法 (teaching-learning-based optimization, TLBO) 是 Rao 等 人 种 算法 进行 了 混合 ; Yong 等 人 外 将 TLBO 算法 作 》 
5] 受 课堂 教学 过 程 启 发 在 2011 年 提出 的 一 种 新 的 群体 智能 优 和 声 的 一 种 策略 ， 将 两 种 算法 进行 了 混合 ， 并 应 用 于 非 负 线性 
化 算法 。TLBO 算法 很 强 的 收敛 能 力 ， 己 经 被 证 明 在 很 多 最 小 二 乘 问题 的 求解 ; Ouyang 等 人 tM 使 用 HS 的 即兴 创作 过 程 
标准 优化 函数 及 工程 问题 上 相 较 于 其 他 许多 全 局 优化 算法 具有 来 优化 “ 教 阶段 ”的 最 优 个 体 ， 而 提出 了 HTLBO-HS 算法 。 现 
更 好 的 表现 ， 是 一 种 有 效 的 全 局 优化 方法 。 但 研究 揭示 [9 该 算 有 这 些 混 合算 法 在 一 定 程度 上 提高 了 算法 的 性 能 ， 但 并 没有 从 
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录用 定稿 


根本 上 交 服 它们 的 固有 缺点 。 使 用 合适 的 融合 方式 将 两 种 算法 
进行 混合 ， 才 能 充分 发 挥 它们 各 自 的 优点 ， 取 得 优势 互补 的 效 
果 。 为 此 ， 本 文 提 出 了 一 种 和 声 教 与 学 混合 算法 (hybrid 
algorithm based on HS and TLBO , HHSTL). HHSTL 算法 依据 
种 群 活跃 率 及 种 群 最 优 个 体 更 新 率 这 两 种 机 制 动 态 地 确定 HS 
算法 或 TLBO 算法 在 下 一 周期 中 作为 种 群 更 新 方式 的 比例 。 利 
用 国际 上 通用 的 16 个 不 同类 型 的 Benchmark 测试 函数 0521 进 
行 仿真 检验 , 并 与 七 种 优化 算法 (DIHS03 IHSU4, RCGHSU!, 
TLBO, ITLBO!9, MATLBO!", TLBO-GCU*9) 进行 了 比较 ， 
检验 了 所 提 算 法 的 性 能 。 
1 ”和 声 搜索 算法 
1.1 标准 和 声 搜索 算法 

标准 和 声 算法 的 基本 步骤 中 如 下 : 

a) 参数 设置 。 

设置 问题 的 维 数 p ， 和 声 记忆 库 的 大 小 HMS ， 最 大 迭代 
WT Ü ERRER HMCR ， 局 部 微调 概率 PAR, Jani 
微调 步 长 Bw 。 
bo 随机 初始 化 和 声 记 忆 库 HM o 


x x x Xp 
2 2 2 2 
HM 4| * x x; X5 
HMS HMS HMS HMS 
x x x Xp 


c) 创作 一 个 新 的 和 声 。 

每 次 迁 代 可 产生 一 个 新 和 声 xn 。 新 和 声 的 每 一 个 音调 
x (i12, D) 按 以 下 三 种 规则 产生 : 以 概率 HMCR 从 和 
声 记忆 库 HM I iH Doors] 中 选择 ; @ 将 四 中 产生 的 


”以 概率 PAR 进行 微调 ， 方 法 如 下 : 
x" = X" + rand : Bw(i) (1) 

加 以 概率 1- HMCR 在 搜索 空间 内 随机 生成 ， 有 具体 如 下 
x^ = xL, + rand : (xU, — xL.) (2) 


d) 更 新 和 声 记 忆 库 。 

对 步骤 c) 中 产生 的 新 和 声 x" 进行 评估 ， 如 果 优 于 HM 
中 的 最 差 和 声 x", WI] xv Ede x, BUF c)。 重 复 步 
又 c) d)， 直 到 满足 终止 条 件 停止。 
1.2 ”改进 的 和 声 算法 

为 改进 和 声 算法 的 性 能 ， 文 献 [14] 提 出 了 一 种 改进 的 和 j 
算法 (IHS )， 算 法 对 参数 PAR 和 Bw 进行 动态 地 调整 。 参 数 
PAR 随和 途 代 的 增加 而 线性 增 大 ， Bw 呈 指 数 地 递减 。 


min 


PAR(t) = PAR n + ME a xt (3) 
In( B's ) 
BW max (4) 


Bw(t) = Bw nax exp( xt) 


其 中 : Ta 表示 最 大 迭代 次 数 ，/ 为 当前 迭代 次 数 ;， PAR na > 
PAR 分 别 为 最 大 、 最 小 微调 概率 ; Bv... 、Bw 分 别 为 最 大 、 


最 小 微调 步 长 。 
2 ” 教 与 学 优化 算法 


21 标准 的 教 与 学 优化 算法 

标准 TLBO 算法 模拟 了 班级 的 教学 和 学 习 行 为 ， 主 要 包括 
“ 教 ”" 阶 段 和 “学 ”阶段 两 个 阶段 ,“ 教 "阶段 , 老师 向 每 位 学 员 传授 
知识 ， 提 高 班级 的 整体 水 平 ;“ 学 ”阶段 ， 学 员 之 间 相 互 学 习 ， 
取长补短 ， 共 同 进步 。 基 本 步骤 由 如 下 : 

a) 参数 设置 。 
NP 表示 班级 的 人 数 (种 群 的 规模 ), D 表示 学 习 的 科 
的 维 数 ), Th 为 最 大 迭代 次 数 。 

b) 在 解 空间 随机 初始 化 NP 位 学 员 : 

Xi=XL+rand:(XL- XU) (5) 
其 中 : j=1,2,…,NP ; rand 为 (0.D 间 的 随机 数 。 
co)“ 教 "阶段 。 


数 (问题 


Xj, = Xda + difference (6) 


difference = rand (1, D) x (x en N nan ) (7) 


其 中 : X44. Xis 分 别 表 示 第 /位 学 员 x3 在 教学 前 、 后 的 知 
RKP; xe 为 老师 ( 即 适应 度 值 最 优 的 学 员 ); 
T, =round( + rand(0,1)) 表示 教学 因子 ， 取 值 为 1 或 2; 


ean 1 d j ` ML Et ` AT 
xU" -aup2, 为 班级 学 员 知 识 水 平 的 均值 。 
j 


适应 度 值 A(X) 反映 了 学 员 X: 的 学 习 成 绩 , 教学 完成 后 ， 
如 果 SALSA ， 则 更 新 个 体 ， 即 X7 =X ， 否 则 不 更 
新 。 

由 “学 "阶段 。 
;从 班级 中 随机 
r 关 7， 进行 交流 学 习 。 
xi d f(X')< fO) 

X}ı +rand:(X!'-X') f(X’)< f(X”) 

学 习 完 成 后 ， 如 果 SALSA ， 则 更 新 个 体 ， 即 
X/-XL,, WA SG. 

e) 如果 满 足 终止 条 件 ， 则 优化 结束 并 输出 最 优 解 ， 否 则 转 
至 c) 继续 。 
22 ” 教 与 学 优化 算法 的 改进 

TLBO 算法 具有 较 强 的 收敛 能 力 ， 但 易于 陷入 局 部 最 优 ， 
尤其 对 中 心 旋转 偏 移 的 复杂 问题 ， 存 在 固有 的 缺陷 。 为 增强 算 
法 的 寻 优 性 能 ， 本 文 做 了 如 下 改进 : 
a) 班级 中 的 最 差 成 绩 往往 容易 拉 低 班级 的 平均 成 绩 , 故 用 
班级 中 的 成 绩 最 差 者 XX" 去 葵 换 班级 平均 水 平 X"…” ， 可 以 
尽快 提升 班级 的 整体 水 平 ， 即 


取 一 个 学 员 xX”，r=142,…,NP， 


AME 
PIN 


(8) 


N 


difference = rand (1, D) x (x sacher TX we) (9) 


b) 学 生 学 习 时 , 不仅 会 在 同学 之 间 进 行 交流 ， 
求教 ， 故 对 学 习 方 式 做 如 下 改进 : 


也 会 向 老师 


201810.00033v1 
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H 的 值 随 欠 代 的 增加 而 逐渐 减 小 ,表明 在 和 迭代 的 前 期 ， 同 
学 间 的 交流 较 多 ， 但 随 着 迭代 的 进行 ， 学 生 会 更 多 地 向 老师 请 


(10) 


(11) 


教 。 王 表示 学 习 因 子 ， 取 值 于 (0,1) 间 。 
c) 学 生 在 学 习 之 余 ， 也 会 根据 自己 当前 的 学 习 状 况 ， 挑 选 
部 分 感 兴趣 的 课程 进行 “自学 ”， 从 而 提升 自己 的 学 业 水 平 。 具 


体 如 下 : 对 第 i 门 课程 , 产生 一 个 (0.D 间 的 随机 数 rand , 用 于 
不 同 的 “自学 ”方式 。 


Xait2x(rand—0.5):4 if rand « SCR 
XD =) X +randx (x -x ) elseif rand < SDR (12) 
Xx otherwise 
P 2 
WEIT T 
max ZA 


i=l,..D,j=1,...,NP , 其中: SCR (self-culture rate, SCR) 


为 自学 概率 ; SDR Cself-developmentrate, SDR) 为 自我 拓展 


Ha. 表示 自学 步 长 ， 和。 、 


分 别 为 最 大 、 最 小 步 长 。 


对 每 一 门 课程 , 学 生 以 SCR 的 概率 进行 自我 调整 , 以 SDR 
的 概率 在 解 空间 随机 搜索 ， 进 行 自我 拓展 学 习 。 为 平衡 算法 在 
搜索 过 程 中 的 空间 扰动 能 力 和 局 部 寻 优 能 力 ， 设 计 自 学 步 长 随 
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3 “和 声 ”-' 教 与 学 "混合 优化 算法 


对 不 同 的 优化 目标 ，HS 和 TLBO 算法 所 表现 的 优化 性 能 
各 不 相同 ， 为 充分 发 挥 它们 各 自 的 优点 ， 更 有 效 地 解决 不 同类 
型 的 优化 问题 ， 本 文 设计 了 一 种 动态 选择 策略 ， 利 用 种 群 活 跃 
率 (AR) 及 种 群 最 优 个 体 更 新 率 (BR) 两 种 机 制 对 两 种 算法 进行 
了 融合 ， 动 态 地 确定 在 下 一 周期 中 选择 HS 算法 或 TLBO 算法 
作为 种 群 更 新 方式 的 比例 。 
3.1 动态 选择 策略 

算法 每 经 过 T〈 选 择 周期 ) 次 迭代 后 使 用 DP《〈 动 态 选 择 概 
AE) 来 重新 确定 两 种 算法 的 选择 比例 。 在 第 k 个 选择 周期 ， 每 
次 迭代 首先 生成 一 个 随机 数 rs(0.D) 。 如果 re < pp , 选择 
HS 算法 来 更 新 种 群 ， 否 则 ， 选 择 TLBO 算法 来 更 新 种 群 。 

种 群 活跃 率 CAR) = 更 新 成 功 个 体 数 /产生 新 个 体 总 数 

种 群 最 优 个 体 更 新 率 (BR)= 种 群 最 优 个 体 更 新 次 数 /所 产 
生 种 群 的 数目 


i 


T 
Mks, App" eos, Stones [e 


, +, GURIA 
歼 延 海 ， 等 : 动态 选择 策略 的 和 声 教 与 学 混合 算法 
AR? _ c(k) “° _ c,(k) 
H  NPxTxDP' ' | NPXTx(ü- DP*)* 
(k) b (k) BR” E b,(k) 
" TxDpes 7o Txü-Dbpeuy 
ARO BR 
(k) _ 
odd ARD + ARP BRI + BRO 09 
H T H T 


其 中 : Ta 为 最 大 迭代 次 数 ，ci(k) 、c,(k) 表示 第 k 个 选择 周期 
HS 和 TLBO 算法 更 新 成 功 的 个 体 数 


内 分 别 使 用 HS 和 TLBO 算法 更 新 


内 分 别 使 用 


b,(k) 表示 第 大 个 选择 周 其 


种 群 最 优 个 体 的 次 数 ; 


Ser 


表示 第 大 个 周期 内 使 


J HS 和 TLBO 5i Hr y 4 


AR, ~ 


WESE: DPO 为 第 k 个 周 


在 第 一 个 周期 ， 
在 后 再 
比例 , 


i 的 每 个 周期 ， 使 月 
那 种 算法 在 迭代 过 程 中 更 活跃 ， 


体 的 数 


越 大 ， 被 选择 的 概率 也 就 
越 多 ， 表 明 算 法 越 活跃 ， 


BU. XE 


法 活跃 并 不 


一 代 种 群 的 更 新 方式 ， 


定 能 使 得 种 群 
群 活跃 率 和 最 优 个 体 更 新 率 共同 确定 的 动态 选择 概率 来 决定 新 
也 确定 两 种 算法 的 混合 比例 以 适应 


的 最 优 


动态 


AR 分 别 表示 第 大 个 周 


E 种 群 的 活跃 率 ， BR 、 


7i 


3 b (k) ` 


期 内 


2 


BR, 分别 


] HS 和 TLBO 算法 所 产生 种 群 的 最 优 个 
期 的 动态 选择 概率 。 

种 算法 以 相同 的 概率 被 选择 , 公平 竞争 ; 
动态 选择 概率 来 决定 两 种 算法 的 选择 
上 且 对 算法 寻 优 的 贡献 


群 更 新 成 功 个 


适用 于 所 求解 的 问题 ， 但 算 
解 得 到 更 新 ， 


所 以 


由 种 


E 


不 同 的 优化 对 象 ， 能 够 充分 发 挥 两 种 算法 的 长 处 ， 提 高 算法 的 


寻 优 能 力 。 针 对 不 同 的 优化 
策略 能 够 根据 两 种 算法 的 


HERI) 


问题 和 同一 问题 的 不 同 阶段 ， 选 择 
态 地 进行 调节 ， 充 分 发 挥 


HS 算法 优秀 的 全 局 搜索 能 力 和 TLBO 算法 的 快速 收敛 能 


3.2 ”算法 流程 


HHSTLBO 算法 的 


L 体 流程 如 图 1 所 示 。 


ETD 
v 
初始 化 参数 和 种 


Y 


tz LDPP-0.5,s-0, k=0 


n pees = 
ZÉ —— rand«pp9 一 二 
Le 


Y Y 
使 用 IHS 算 法 NP 次 ， 使 用 改进 的 TLBO 算 法 一 次 ， 
更 新 a1(k) 和 bi(k) 更 新 aa(k) 和 bs(k) 
— S=S+1 四 —— 
是 X. 
上 S<T pe 
y 8 
s=0,k=k+1 


— 


Y 
按 式 (14) 计算 DP = 


图 1 HHSTL 算法 流程 
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4 ”实验 仿真 测试 


为 检验 HHSTL 算法 的 性 能 ， 将 其 与 另外 七 种 优秀 算法 
(DIHS, HIS. RCGHS. TLBO. ITLBO, MTLBO. TLBO-GC) 
在 16 个 不 同类 型 的 标准 测试 函数 上 进行 仿真 测试 。 其 中 ，F1~ 
Fs 是 单 峰 函数 ，F4 是 混合 函数 ，Fs ~ Fo 是 复杂 的 多 峰 函 数 ， 
Fi3~Fi6 是 对 4 个 经 典 的 多 峰 函 数 进行 了 Shift 操作 ， 这 些 类 型 
各 异 的 标准 函数 能 很 好 地 测试 算法 的 通用 性 。 

Sphere Shift Function 


Fi: Sphere Function F2: 
F3 : Zakh Function F4: Bohachevsky Shift Function 
Fs: Ackley Function Fe: Griewank Function 
F7: Levy Function Fs: Schwefel2.22 Function 


Fo: Rastrigin Function Fio: Michalewics Function 
Fu: Schwefel Function Fi2: 
Fi3: Ackley Shift Function 


Fi4: Griewank Shift 


Schwefel 2.26 Function 


Function 

Fıs: Rastrgin Shift Function 

Fie: Schaffer Shift 
41 实验 环境 与 参数 设置 

本 文 所 有 测试 的 环境 为 : 戴尔 PC 机 Intel Core TM i7-4790 
3.6 GHz CPU, 8 GB 内 存 ; Window7 操作 系统 , MATLAB R2014b 
软件 。 为 保证 测试 的 公平 ， 各 种 算法 采用 相同 的 适应 度 函 数 评 
MRZE (5000x 维 数 )。 比 较 算法 的 参数 设置 参考 于 原文 献 ， 设 
具体 如 表 1 所 示 。 


Function 


表 1 参数 设置 


算法 种 群 大 小 其 他 
HMCR = 0.99, PAR ax = 0.99, PAR ,, = 0.15 
DIHS 10 Bw = (xU — xL) / 20, Bw in = (xU — xL) / (le +15) 
S. 2200/ D, Spa =20/ D, fw y = (xU — xL) / (le+4) 
HMCR = 0.95, PAR,,, = 0.99, PAR,,, = 0.45 
IHS 5 
Bw ax = (XU — xL) / 20, Bw „in = 0.0001 
RCGHS 5 / 
TLBO 10 / 
ITLBO 40 / 
MTLBO 10 SR, 7 0.8, SR, = 0.3, SFR = 0.35 
TLBO-GC 10 / 
HMCR = 0.99, PAR „ = 0.99, PAR,, =0.1 , 
下 =035，Bw = A = GU — xL) /10 
HHSTL 10 


Bw, = Anin = (xU — xL) / (le +15) 


SCR = 0.3,SDR = 0.03,T =150 


42 实验 结果 与 分 析 

表 2 统计 了 16 个 测试 函数 在 D=100 时 的 30 次 独立 实验 结 
R, 包括 每 个 算法 的 目标 函数 适应 度 值 的 最 优 值 (Best)、 平均 值 
(Mean)、 最 差 值 (Worst) 和 方差 (Std)。 所 测试 的 函数 中 除 Fio 


以 外 ,理想 最 优 值 均 为 0， 所 以 均值 越 接 近 0, 表明 算法 收敛 精 
度 越 高 ， 而 方差 越 小 ， 表 明 运 算 结果 的 浮动 范围 越 小 ， 算 法 越 
从 表 2 可 以 看 出 ，RCGHS、TLBO、ITLBO 及 TLBO-GC 
算法 在 四 个 函数 Fi、F2、Fe、Fs 能 达到 理论 最 优 值 , 但 在 Fi1~Fi6 
上 效果 明显 很 差 ， DIHS 算法 在 七 个 函数 F4、 Fs. Fe. Fo. Fins 
Fi, Fia 能 取得 最 优 值 ，IHS 算法 在 三 个 函数 Fio、Fun、F1lz 能 取 
GRRE, BENER Fi Fa Fo Fs 上 收敛 精度 都 不 高 ; 
而 本 文 所 提 的 HHSTL 算法 在 16 个 测试 函数 上 的 平均 值 和 方差 
均 能 获得 最 优 值 ， 说 明 所 提 算 法 综合 了 HS 和 TLBO 两 种 算法 
的 优点 , 具有 较 强 的 优化 能 力 。 为 进一步 分 析 实 验 结果 , 图 2、 
3 给 出 了 八 种 算法 30 次 独立 运行 的 平均 收敛 曲线 图 和 最 优 值 分 
HAR (D=100 维 )。 从 收敛 曲线 图 不 难看 出 ， 本 文 算法 对 两 个 
复杂 优化 函数 (Rastrigin Shift Function 、Schaffer Shift 
Function) 的 收敛 曲线 呈 下 降 趋势 , 收敛 效果 明显 优 于 其 他 算法 ， 
说 明 本 文 算法 的 全 局 搜索 能 力 较 强 ,不易 陷 入 局 部 最 优 。 从 统 
计 盒 图 可 以 看 出 ， 本 文 算 法 的 30 次 最 优 值 的 分 布 基本 呈 一 条 
直线 ， 说 明 本 文 算法 具有 很 好 的 稳定 性 
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图 3 最 优 值 分 布 统计 盒 图 

为 检验 本 文 算法 与 比较 算法 之 间 的 差异 性 ， 表 3 列 出 了 本 
文 算 法 与 七 种 比较 算法 在 16 个 测试 问题 上 30 次 独立 运行 结果 
的 Wilcoxon 符号 秩 检验 的 P-value fj, *—". “+”, “~” 4 5l 
表示 比较 算法 与 本 文 算法 的 成 绩 相 比 是 差 、 好 和 相当 (NaN K 
示 算 法 成 绩 无 差异 )。 

从 表 3 可 以 看 出 ， 七 种 比较 算法 在 测试 问题 上 没有 优 于 本 
文 算法 的 , DIHS 在 七 个 问题 上 与 本 文 算法 相当 , IHS 算法 在 六 
个 问题 上 与 本 文 算法 相当 ，RCGHS、TLBO 和 ITLBO 算法 均 
在 五 个 问题 上 与 本 文 算法 相当 。 由 此 可 见 ， 本 文 算 法 与 比较 算 
法 在 大 多 数 测试 问题 上 存在 显著 性 有 明显 的 统计 意义 。 
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录用 定稿 RE, 择 策 学 混合 算法 
表 2 算法 30 次 独立 运行 结果 比较 (D=100) 
ALG Best Mean Worst Std ALG Best Mean Worst Std 
DIHS 2.88E-20 3.15E-20 3.61E-20 2.22E-21 DIHS 1.03E-17 1.18E-17 1.30E-17 7.73E-19 
IHS 9.26E-27 1.03E-26 1.12bE-26 6.56E-28 IHS 3.44E-24 4.00E-24 4.77E-24 4.47E-25 
RCGHS 0 0 0 0 RCGHS 0 0 0 0 
TLBO 0 0 0 0 TLBO 0 0 0 0 
B ITLBO 0 0 0 0 ITLBO 0 0 0 0 
MTLBO 2.24E-84 2.62bE-82 1.86E-81 5.75E-82 MTLBO 1.07E-81 1.48E-80 4.72E-80 1.46E-80 
TLBO-GC 0 0 0 0 TLBO-GC 0 0 0 0 
HHSTL 0 0 0 0 HHSTL 0 0 0 0 
DIHS 2.41E+02 2.80E+02 3.18E+02 3.12E+01 DIHS 0 0 0 0 
IHS 1.40E+02 1.65E+02 1.98E+02 1.94E+01 IHS 0 3.26E-04 3.26E-03 1.03E-03 
RCGHS 3.57E-13 3.87E+14 7.43E+14 2.33E+14 RCGHS 4.87E+03 4.92E+03 4.99E+03 3.50E+01 
TLBO 1.40E-254 9.32E-235 9.31E-234 0 TLBO 5.27E+01 1.62E+02 2.64E+02 7.58E+01 
is ITLBO 1.73E-110 4.17E-102 4.15E-101 1.31E-101 m ITLBO 2.41E+01 3.30E+01 4.29E+01 5.94E+00 
MTLBO 1.31E+02 1.91E+02 2.57E--02 3.41E+01 MTLBO 0 1.80E-14 1.28E-13 3.98E-14 
TLBO-GC 1.63E-179 8.51E-153 8.50E-152 2.69E-152 TLBO-GC 2.98E+03 3.40E+03 3.69E+03 2.28E+02 
HHSTL 0 1.49E-262 1.49E-261 0 HHSTL 0 0 0 0 
DIHS 4.33E-10 4.48E-10 4.67E-10 1.10E-11 DIHS 0 0 0 0 
IHS 2.28E-13 2.5]E-13 2.71E-13 1.15E-14 IHS 0 7.40E-04 7.40E-03 2.34E-03 
RCGHS 4.44E-15 4.80E-15 7.99E-15 1.12E-15 RCGHS 0 0 0 0 
TLBO 4.44E-15 6.93E-15 7.99E-15 1.72E-15 TLBO 0 0 0 0 
FS ITLBO 4.44E-15 6.93E-15 7.99E-15 1.72bE-15 ITLBO 0 0 0 0 
MTLBO 5.06E-14 5.84E-14 6.48E-14 4.37E-15 MTLBO 0 4.395E-13 4.394E-12 1.39E-12 
TLBO-GC 4.44E-15 6.57E-15 7.99E-15 1.83E-15 TLBO-GC 0 0 0 0 
HHSTL 8.88E-16 1.24E-15 4.44E-15 1.12E-15 HHSTL 0 0 0 0 
DIHS 5.29E-20 6.02E-20 6.44E-20 3.77E-21 DIHS 1.085E-09 1.135E-09 1.178E-09 2.927E-11 
IHS 1.68E-26 2.01E-26 2.38E-26 2.45E-27 IHS 7.359E-13 7.83E-13 8.492E-13 3.617E-14 
RCGHS 9.39E--00 9.49E--00 9.57E+00 5.01E-02 RCGHS 0 0 0 0 
TLBO 4.15E+00 5.18E+00 6.18E+00 5.52E-01 TLBO 0 0 0 0 
à ITLBO 4.12E+00 4.62E+00 5.55E+00 4.47E-01 ITLBO 0 0 0 0 
MTLBO 1.24E-23 5.67E-22 3.11E-21 9.24E-22 MTLBO 4.207E-46 3.448E-45 1.189E-44 4.129E-45 
TLBO-GC 7.92E+00 8.21E+00 8.55E+00 2.26E-01 TLBO-GC 0 0 0 0 
HHSTL 4.86E-27 7.15E-27 9.04E-27 1.42E-27 HHSTL 0 0 0 0 
DIHS 0 0 0 0 DIHS -9.93E-*01 -9.93E+01 -9.91E+01 5.48E-02 
IHS 0 3.04E-03 1.30E-02 5.05E-03 IHS -9.94E+01 -9.94E+01 -9.93E+01 4.35E-02 
RCGHS 0 0 0 0 RCGHS -4.09E+01 -3.80E+01 -3.06E+01 3.05E+00 
TLBO 0 8.56E+00 8.56E+01 2.71E+01 TLBO -9.03E+01 -7.90E+01 -5.07E+01 ].12E*01 
"m ITLBO 0 0 0 0 SM ITLBO -8.48E+01 -7.83E+01 -6.64E+01 6.04E+00 
MTLBO 4.80E+00 8.77E+00 1.20E+01 1.97E+00 MTLBO -9.65E+01 -9.61E+01 -9.57E+01 2.71E-01 
TLBO-GC 0 ].12E*02 5.00E--02 1.88E+02 TLBO-GC -4.79E+01 -4.27E+01 -3.70E+01 3.33E+00 
HHSTL 0 0 0 0 HHSTL -9.94E+01 -9.94E+01 -9.93E+01 4.24E-02 
DIHS 1.31E-10 1.31E-10 1.31E-10 0 DIHS 1.31E-10 1.31E-10 1.31E-10 0 
IHS 1.31E-10 1.31E-10 1.31E-10 0 IHS 1.31E-10 1.31E-10 1.31E-10 0 
RCGHS 3.74E+04 3.84E+04 3.94E+04 7.30E+02 RCGHS 3.72E+04 3.82E+04 3.90E+04 6.13E+02 
TLBO 1.52E+04 1.95E+04 2.28E+04 2.52E+03 TLBO 1.09E+04 1.89E+04 2.55E+04 4.90E+03 
EN ITLBO 1.77E+04 1.94E+04 2.19E+04 1.19E+03 ITLBO 1.41E+04 1.74E+04 2.04E+04 1.66E+03 
MTLBO 1.60E-10 2.68E-10 4.58E-10 9.23E-11 MTLBO 1.60E-10 2.44E-10 4.07E-10 7.87E-11 
TLBO-GC 2.65E+04 2.84E+04 3.20E+04 1.58E+03 TLBO-GC 2.55E+04 2.74E+04 2.91E+04 1.19E+03 
HHSTL 1.31E-10 1.31E-10 1.31E-10 0 HHSTL 1.31E-10 1.31E-10 1.31E-10 0 
DIHS 4.19E-10 4.44E-10 4.61E-10 1.38E-11 DIHS 0 0 0 0 
IHS 2.42E-13 2.54E-13 2.70E-13 1.01E-14 IHS 0 I.11E-17 1.11E-16 3.51E-17 
RCGHS 2.06E+01 2.08E+01 2.10E+01 9.15E-02 RCGHS 2.75E+03 3.19E+03 4.68E+03 7.59E+02 
TLBO 1.89E+01 1.91E+01 1.93E+01 1.33E-01 TLBO 3.26E+00 6.89E+01 1.45E+02 4.95E+01 
ni ITLBO 1.77E*01 1.84E+01 1.90E+01 4.06E-01 idi ITLBO 1.03E-02 8.68E-01 4.51E+00 1.44E+00 
MTLBO 3.13E-11 1.13E-10 2.36E-10 6.76E-11 MTLBO 4.44E-16 7.58E-04 7.41E-03 2.34E-03 
TLBO-GC 1.99E+01 2.01E+01 2.04E+01 1.82E-01 TLBO-GC 1.62E+03 1.90E+03 2.12E+03 1.60E+02 
HHSTL 1.99E-13 2.18E-13 2.45E-13 1.01E-14 HHSTL 0 0 0 0 
DIHS 0 3.33E-03 3.19E-02 1.01E-02 DIHS 2.38E-03 6.59E-03 1.77E-02 6.73E-03 
IHS 0 8.42E-03 5.75E-02 1.81E-02 IHS 4.05E--00 9.30E--00 1.84E*01 5.41E+00 
RCGHS 1.77E+03 1.84E+03 1.98E+03 6.08E+01 RCGHS 8.09E+02 8.21E+02 8.48E+02 1.24E+01 
TLBO 7.64E+02 8.73E+02 9.64E+02 6.52E+01 TLBO 5.90E+02 6.12E+02 6.34E+02 1.72E+01 
m ITLBO 6.69E--02 7.52E*02 8.14E*02 5.03E-*01 is ITLBO 5.12E*02 5.42E*02 5.71E*02 2.01E-*01 
MTLBO 2.00E--00 5.59E--00 1.10E+01 2.75E+00 MTLBO 8.60E+01 1.35E+02 1.80E+02 3.96E+01 
TLBO-GC 1.25E+03 1.32E+03 1.45E+03 6.40E+01 TLBO-GC 6.38E+02 6.92E+02 7.38E+02 2.75E+01 
HHSTL 0 5.86E-15 3.91E-14 1.32E-14 HHSTL 4.59E-05 3.32E-03 1.67E-02 5.33E-03 
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录用 定稿 RE, F: 动态 选择 策略 的 和 声 教 与 
表 3 Wilcoxon 符号 秩 检验 的 P-value (D-100) 4.3 ”选择 周期 T 的 分 析 
函数 DIHS IHS RCGHS TLBO  ITLBO MTLBO TLBO-GC HERAA T RAR SERRA SIS PE 02 E T fü 
Fl 639E-05 6.39E-05 NaN NaN NaN  639E-05  NaN 定 下 一 周期 中 两 种 算法 作为 种 群 更 新 方式 的 比例 。 为 检验 选择 
F2 6.39B-05 6.39E-05 NaN NaN NaN 6.39E-05  NaN 周期 工 对 算法 的 影响 , 在 六 个 测试 函数 上 对 工 的 不 同 取 值 进行 
F3 1.79E-04 1.79E-04 1.79E-04 1.79E-04 1.79E-04 1.79E-04 1.79E-04 
F4 6.16E-05 1.71E-03 1.78E-04 1.78E-04 1.78E-04 2.72E-01 1.78E-04 二 
F5 8.74E-05 8.45E-05 1.24E-04 9.14E-05 9.14E-05 7.57E-05 1.17E-04 应 度 值 的 乎 均值， 图 4 给 出 了 不 同 工 值 下 30 次 测试 的 分 布 盒 
F6 NaN — 3.68E-0] NaN NaN NaN 7.79E-02 NaN 图 。 
F7 1.83E-04 1.83E-04 1.83E-04 1.83E-04 1.83E-04 1.83E-04 1.83E-04 从 表 4 可 以 看 出 ， 除 函数 Griewank shift Function 4h, 多数 
F8 6.39E-05 6.39E-05 NaN NaN NaN 。 6.39E-05 NaN 测试 函数 在 不 同 T 值 下 得 到 的 平均 适应 度 值 差 异 不 是 很 大 ， 说 
F9 NaN 1.49E-02 NaN 3.68E-01 NaN 6.39E-05 7.79E-02 明 所 提 算 法 对 此 参数 不 大 敏感 。 根 据 表 4 和 图 4 的 结果 ， 综 合 
Fl0 1.71E-03 2.12E-01 1.83E-04 1.83E-04 1.83E-04 1.83E-04 1.83E-04 考虑 解 的 精度 及 分 布 统计 情况 ， 本 文 取 T=150， 
FH NaN NaN 6.39E-05 6.39E-05 6.39E-05 6.29E-05 6.39E-05 
F12 NaN NaN 6.39E-05 6.39E-05 6.39E-05 6.34E-05 6.39E-05 — ] C 
Fl3 — L71E-03 2.12E-01 1.83E-04 1.83E-04 1.83E-04 1.83E-04 L83E-04 E | | o 
F14 NaN 3.68E-01 6.39E-05 6.39E-05 6.39E-05 6.39E-05 6.39E-05 了 H |] | _ | SDC | afl- 
F15 4.56E-01 837E-03 1.11E-04 1.11E-04 1.11E-04 1.11E-04 1.11E-04 i . |] | | I | | : | | u | | |! 
Fl6 121E-01 L83E-04 1.83E-04 1.83E-04 L83E-04 1.83E-04 L83E-04 | a f n | | | Í H | Hn Hi 
E 9 10 "N 11 11 14 13 [H ELA " H R Í | | | 
+ 0 0 0 0 0 0 0 C NCTCETOETOETTNCT TNT RT RT 
~ 7 6 5 5 5 2 3 
图 4 不 同 T 值 下 的 最 优 值 分 布 统计 盒 图 
表 4 不 同 T 值 下 30 次 独立 实验 的 最 优 目 标 函 数 的 平均 适应 值 
Function T=5 T-50 T-100 T-150 T-200 T-300 T-500 T-1000 
Schwefel Function 8.37E-12 728E-12 728E-12 728E-12 7.64E-12 8.37E-12 7.64E-12 7.64E-12 
Levy Function 1.07E-29 9.05E-30 1.21E-29 1.57E-30 1.63E-29 3.06E-29 1.75E-29 1.32E-29 
Schwefel2.26 Function 8.37E-12 728E-12 728E-12 728E-12 7.64E-12 8.37E-12 7.64E-12 7.64E-12 
Ackley shift Function 7.18E-14 6.59E-14 7.59E-14 6.25E-14 6.79E-14 6.32E-14 6.84E-14 6.ATE-14 
Rastrigin shift Function 3.55E-16 2.49E-15 2.49E-15 3.55E-16 2.84E-15 2.13E-15 3.55E-16 1.07E-15 
Griewank shift function 0 4.93E-03 0 0 6.40E-03 8.37E-03 3.94E-03 7.87E-03 
4.4 动态 选择 策略 分 析 SchwefeD.26 Function 这 样 的 复杂 多 峰 函 数 ， 由 于 存在 多 个 局 
为 进一步 分 析 本 文 算法 所 设计 的 动态 选择 策略 ， 分 别 选 择 ” 部 极 小 值 点 , 进化 的 前 期 需要 进行 全 局 搜索 , HS 算法 因为 较 强 
了 单 峰 和 复杂 多 峰 函 数 进行 了 测试 。 图 5 给 出 了 两 种 不 同类 型 ” 的 全 局 搜索 能 力 而 获得 较 高 的 选择 比例 ， 能 快速 定位 最 优 解 的 
函数 的 动态 选择 概率 变化 曲线 图 。 有 效 区 域 。 随 着 和 迭代 的 进行 ， 两 种 算法 的 选择 概率 慢 慢 趋 近 ， 
ETI ， — ————— 交替 变化 ， 在 迭代 的 后 期 ，TLBO 算法 的 选择 概率 增 大 ， 主 要 
"| Va Nw, 3M — e| 来 进行 精细 搜索 ， 提 高 解 的 精度 。 选 代 过 程 中 ， 两 种 算法 都 保 
ol — e AL 以 a" 持 着 一 定 的 选择 比例 ，HS 算法 能 够 维持 种 群 的 多 样 性 ， 而 
| š jJ | dr 一 一 一 TLBO 算法 能 够 加 快 收敛 速度 ， 提 高 解 的 精度 。 
jh hi "yd Vra A al Mugs ^ | | | 5 RE 
idi LE 为 提高 对 多 种 不 同类 型 问题 的 优化 能 力 ， 本 文 提出 了 一 种 
图 5 动态 选择 概率 变化 曲线 图 和 声 教 与 学 混合 优化 算法 一 一 HHSTL。 算法 使 用 由 种 群 活跃 率 
从 图 5 中 可 以 看 出 ， 动 态 选择 策略 能 够 根据 优化 问题 的 特 。” ”及 种 群 最 优 个 体 更 新 率 构成 的 动态 选择 概率 将 两 种 算法 进行 了 
点 和 迭代 过 程 中 两 种 优化 算法 的 具体 表现 动态 地 选择 合适 的 更 。 有 效 地 融合 。 通 过 在 16 个 不 同类 型 的 Benchmark 函数 上 对 
新 策略 来 满足 不 同 进化 阶段 的 需求 。 对 Sphere Function 这 样 的 HHSTL 算法 及 其 他 七 种 优秀 算法 进行 仿真 比较 , 统计 结果 显示 ， 
单 峰 函数 ，TLBO 算法 因为 良好 的 收敛 能 力 ， 因 而 在 整个 优化 ”所 提 算 法 汲取 了 HS 和 TLBO 两 种 算法 的 长 处 ， 优 势 互补 ， 求 
过 程 中 占据 主导 地 位 ， 选 择 概率 明显 高 于 HS 算法 ; 对 解 精度 高 ， 稳 定性 好 ， 能 在 较 多 不 同类 型 的 优化 问题 上 取得 良 


好 的 效果 。 但 因为 HHSTL 算法 是 由 IHS 4 TLBO 两 种 算法 混 
合 而 成 ， 参 数 较 多 ， 需 要 进行 大 量 的 测试 。 在 后 续 研 究 中 ， 可 
以 设计 有 效 的 参数 自 适应 策略 来 进行 改进 ， 也 可 以 将 该 算法 应 
日 到 实际 的 工程 优化 问题 中 ， 检 验 其 具体 实用 性 。 
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